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ABSTRACT 

The formation of massive stars exceeding 10 M Q usually results in large-scale molecular outflows. 
Numerical simulations, including ionization, of the formation of such stars show evidence for ionization- 
driven molecular outflows. We here examine whether the outflows seen in these models reproduce 
the observations. We compute synthetic ALMA and CARMA maps of CO emission lines of the 
outflows, and compare their signatures to existing single-dish and interferometric data. We find that 
the ionization-driven models can only reproduce weak outflows around high-mass star-forming regions. 
We argue that expanding H n regions probably do not represent the dominant mechanism for driving 
observed outflows. We suggest instead that observed outflows are driven by the collective action of 
the outflows from the many lower-mass stars that inevitably form around young massive stars in a 
cluster. 



1. INTRODUCTION 

Jets and outflows always accompany the star forma- 
tion pro cess, both in the low-mass and in the high-mass 
regime dCabrit k Andre||1991[ |Bachiller||1996l |Reipurth| 
& Bal lyjiWmt |Shepherd|2l)U5l |Beut,her k Shepherd|iW()51 
Arce et al.||2007| |Bally et al.||2007[ |Bally||20071 |2008p . 
Though outflows from low-mass stars are m general 
understood as resulting from magnetically-driven disk 
winds during the accretion process, the formation of 
high-mass outflows from ionizing stars with masses in 
excess of lOAf is less understood, as we argue in the 
next Section, both from the obse rvati onal (Section 2.1 1 
as well as the theoretic al (Section 2.2 ) point of view 



The recent finding by Peters et aT " P010a[ ) that ioniza- 
tion feedback from massive protostars can drive bipo- 
lar molecular outflows has suggested that these could 
explain the observations of high-mass outflows. In the 
present paper, we compare the morphology of these 
molecular outflows with synthetic CO J = 2 — 1 obser- 
vations (Section [3]). We demonstrate that the synthetic 
observations fail to consistently reproduce observations 
of molecular outflows in Section [Ij In Section [5j we dis- 
cuss the implications and limitations of this result, and 
speculate on the actual origin of the observed outflows. 

2. OBSERVATIONS AND THEORY OF HIGH-MASS 
OUTFLOWS 

2.1. Observations of High- Mass Outflows 
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Bipolar outflows from high-mass stars are, in many 
ways, scaled-up versions of their lower-mass counter- 
parts. They do, however, carry significantly mor e mo- 
mentum (e.g. |Beuther et "al]|2002| [Wu et al.||2004j) and 
are often less collimated (e.g. peuther fe Shepnerd| 2005 ) . 
Outflows are launched from accretion disks around 



young stellar objects (e.g. |Ray et aT 20071. Disk-like 



mass 


protostars (Chini et al. 


|2004l |2006[ iBelran et al. 


2004 , 2006j|Jiang et al. 2008; 


D avies et al.||2010| |Zap at a 


et al. 


2010; Motogi et al. 


|201 


ip, and even in situations 



(Howard et al. 


1997; Zhang et al.||1998| |Klaassen et al. 


2009; Ualvan-fv! 


ladrid et al. 2009; 




ruya et al.||2008, 2011) and ionized ( 


Garay et al.||1986; 


Sewilo et aL 


pDDSn^-eto & Klaassen 


20U8|) gas can still 



be detected. In many of these regions, bipolar molecular 
outflows originating from the center of rotation and ori- 
ented perpendicular to the rotation have been detected 
as well, suggesting a causal connection between rotation 
and outflow. 

If high-mass outflows were driven magnetically like 
their lower-mass analogs, then one would expect a corre- 
sponding magnetic field structure along the outflow. In- 
deed, recent observations suggest a correlation between 
magnetic field orientatio n and the presence of outflows 
from hi; 



gh-mass sources dGirart et al.j [20091 |Tang et al- 
ii [20101 rvTemmings et al.||2010| |Beuther et al. 



2009a b 

2010| |i jurcis et al.||201l[ ). However, whether these maser 



and dust polarization measurements indicate that the 
outflows organize the magnetic field or the magnetic field 
collimates the outflows is less clear. 
The momentum and energy imparted to molecular out- 



of the powering source (Richer et al.|2000 


Beuther et al. 


2002 


Arce et al. 2007 


Lopez- Scpulcrc et a 


.|2009|. Once 



protostar becomes luminous enough to 
ionize its surroundings, the interaction between outflow 
and H n region further complicates the picture. Observa- 
tions sug gest that the outflows originate within the ion- 
ized gas ( |De Pree et al.|1994[ |Galvan-Madrid et al.|2009 



Klaassen et al. in prep.), and that they can continue 
into a molecular outflow beyond the ionization bound- 



2 



ary ( e.g. K3-50A; De Pree et al.||1994 and Howard et al 



1997]). 

In order to make statistical arguments about the na- 
ture of outflows from protostars of at least intermediate 
mass, we r ely on single dish surveys to cha racterize their 
properties. Shepherd fc Churchwell| (1996) looked at ten 
sources with bolometric luminosities exceeding 10 2 Lq 
and only found 5 bipolar outflows at a resolution of 60", 
noting that for the other sources, two had multiple veloc- 
ity components, and the other three ha d signal to noise 
ratios too low to draw conclusions from. |L6pez-Sepulcre 
(2009), at a spatial resolution of 11", looked at 



et al. 



higher luminosity sources (£bol > 2x 10 4 L ) and found 
outflows for each of the nine sources they mapped, and 
strong line wings in the other two sources. Mo st of these 
last sources are associated with H II regions. |Qin et aT] 



(2008) studied the gas motions from 15 high-mass star 



forming regions with ultracompact H n regions and found 
evidence for outflows towards 10 of them. However, at 
a spatial resolution > 1', this detection fraction should 
be used as a lower limit. The five sources without out- 
flows were amongst the nearest sources in their survey, 
suggesting that distance did not play a role in the non- 
detections. In each of these studies, outflow masses range 
from 3 to 10 3 Mq, and their energies are of the order 10 46 
erg or greater. 

With interferometers one gains a better understand- 
ing of which source is responsible for the large-scale out- 
flows seen in the single dish surveys. With current in- 
struments we are limited in the number of sources for 
which detailed analyses can be performed. Nonetheless, 
in-depth observations have been made by many groups 
in an attempt to better understand the powering sources 
and their relationships with their environments. Gener- 
ally, these high resolution (< 1") studies show that the 
outflow originates at the brightest contin uum object in 
the field. These sources include: G29.96 (Beltran et al. 



2011 



Klaassc 



|, G24.78A (|Furuya et a l1[2002l W51e2 (IKetofc ; 
ssen|2008l|Shie t al. 20l O|7G3^(|A cord et al.|l997l 



Watson et al.||2007| |Hunter et al.||2008' ), NGC_7538jRS 
1 and 28.2 fKlaassen et al.||201If7G31.41 (|Qlmi et 
1996a|bp , G75.78 and G75.77 ( |Rrrfel fc Liidke||2010| . 



In many of these sources, rotation ot the warm gas sur- 
rounding the H II region has also be en detected perpen- 
dicular to the outflow dire ction (see [B eltran et al. | |2006[ 
2011| Klaasscn et al. 2 0091 |Qlmi eta l. 1996b .^etcTJ" 
'Keto & Klaassenl d2008|) and Klaasscn et al.| (12011 1 



3 poos 

showed that molecular outflows appear at the edges of 
the H II regions already at high velocities, suggesting 
these outflows are powered from within the H II region. 
In fact, the radio recombination line study of De Pree 



et al. (1994) displayed velocity shifts in this ionized gas 



tracer from K3-50A, sugges ting an ionized ou tflow. The 
molecular observations of Howard et al. ( 19971 are consis- 
tent with this outflow penetrating the ionization bound- 
ary and continuing as a molecular outflow. Observa- 
tions of ion ized outflow lobes ha ve also been re ported 



for G48.75 (Johnston et al. 2009). Furthermore, Cesa- 



roni et al. ( 2010p discuss the possibility that the outflow 



in G31.4T o riginates from an ionized jet, and [Guzman 



et al. (2011 1 suggest that an ionized jet in IRAS 16562- 
3959 is the energy source of a molecular outflow in the 
same region. 



2.2. Possible Outflow Driving Mechanisms 

We now consider the problems with existing explana- 
tions of outflow driving. Protostellar outflows around 
low-mass protostars are generally attributed to m agnetic 



fields, either through the launching of disk winds ([B land 



ford fc Payne 



dritz 
netic 



1982 



Pudritz 

Pudritzl 119971 IKrasnopolsky et al.||1999"| IKomgT 
20001 iFendt 



Norman 1983; Ouycd & 
Fu" 



tower flows 



Cemeljic||2002| 
d'l'omisakal 1998 



Fendt 20 06b, mag- 
|2002| |Lynden-Bell 



ell 

14; 



20031 i Matsum oto fc Tom isaka 2004; Mac hida et"aT . 200 
Banerjee fc Pudritz|)2006[ |2007| |Seifried et al.||2012[ ) or 
through the interaction be tween the protostellar m agnc 



tosphere and the disk field flNajita fc Shu|1994"l|Shu et al 
T9941 [20071 |Cai et al.||2008] |Lc 1 



manova et a 



Cai et al. 1 12008 1 |Lovelace et 
■| |2009[ ). Recently, the applicability of 



aijirg 



m iRo- 

oTthese 



ese 

driving mechanism s in the h igh-mass regime has been 
questioned ( Peters et al-[2011 ) on the grounds that severe 
gravitational instability m the accretion flow surrounding 
the massive star destroys the coherent rotational motion 
in the disk plane that is required for all these models 
to work, and the ionization feedback disrupts the mag- 
netic field structure. The statistical distribution of or- 
bital separations wi th a peak at a semimajor axis of 
~ 0.15 AU found by |Kobulnicky fc Fryer] (|2007| ) in the 
massive binary data of Garman y et al.| ( |l980[ ) suggests 
that the inner disk region is gravitationally unstable as 
well, making magnetic jet launching very difficult. Fur- 
thermore, magnetically driven jets around low-mass stars 
may be collima ted at the disk-magnetosphere boundary 
(Lii et al. 2012 ), and massive stars are known to have rel- 
ativ ely weak magn etospheres compared to low-mass stars 
(e.g. Owocki 2009). To fully model this process we need 
to resolve the inner disk region (radii less than 1 AU) 
where magnetic launching is expected. This will require 
high-resolution radiation-magnetohydrodynamical calcu- 
lations that connect the larger-scale gravitationally un- 
stable accretion flow (radii of several thousand AU) to 

t he inner disk region. 

2011[) tried to fill this gap in two- 



Vaidya et al. 



( [20091 

dimensional calculations that include both ideal mag- 
netohydrodynamics as well as the effects of radiative 
pressure in the line-driving approximation. They find 
that initially well collimated outflows can widen up as 
the central star grows in mass and its luminosity in- 
creases. This could explain the observational proposal 
that more evolved outflows around massive stars have 



larger opening an gles than very young ones (Beuther 



& Shepherd 2005). Because of the restricted geometry 
in the Vaidya et al. work, however, these studies can- 
not consider the effects of disk fragmentation, so further 
analysis is required. A more consistent treatment of the 
complex interplay between the fragmentation behavior of 
the inner disk, magne tic fields, and radiative p rocesses 
has been presented by |Commeroon et al. (2011). Their 
three-dimensional radiation-magnetohydrodynamics cal- 
culations indicate that magnetic braking efficiently re- 
moves angular momentum from the disk, leading to high 
accretion rates and luminosities. This heats the disk and 
reduces the degree of fragmentation at small radii. Be- 
cause these calculations are prohibitively expensive, they 
could not follow the evolution into the realm of mas- 
sive stars where ionizing radiation becomes important . 
This regime has been covered by Peters et al. (20111 
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with somewhat lower resolution. They find that mag- 
netic fields cannot comp letely suppress fragmen tation at 
large disk radii (see also Hennebelle et al. 2011 1. 



Additional evidence against magnetic driving mecha- 
nisms comes from the notable lack of observations of ex- 
tremely high- velocity jets associated with high-mass out- 
flows. Since magnetically launched jets form deep within 
the gravitational potential well around the star, and 
the potential well around high-mass stars is deeper than 
around low-mass stars, jets would be expected with ve- 
locities comparable to those of line-driven stellar winds, 
exceeding 1000 km s -1 . Instead, the observed structures 
never propagate faster than jets from low-mass stars. 
However, even if there were very high- velocity jets, they 
might be immediately slowed by interaction with the ac- 
cretion flow, which is expected to be stronger for high- 
mass than for low-mass stars. 

Radiation pr essure on the accretio n flow could also 
drive outflows. Krumholz et al. (2009) used a gray flux- 
limited diffusion method to simulate this effect. They 
found that the radiation pressure creates bubbles around 
the massive star that get destroyed by Rayleigh-Tayor in- 
stabilities. These radiation bubbles are morphologically 
dissimilar to molecular outflows. Additionally, by the 
time the bubbles form, the star would be massive enough 
to ionize the interior of the bubble, so that the bubble 
would become an H II reg ion and not be molecular any- 
more. iKuiper et al . I pOllI ) , using a frequency-dependent 
treatment of the direct stellar radiation combined with 
gray flux-limited diffusion, do not see these Rayleigh- 
Tayor instabilities. As possible reasons for this difference 
they suggest a better modeling of the first absorption of 
photons with their numerical scheme, but it could also 
be that they will see the inst ability once their bubb les 



reach the same size as those of Krumholz et al 



(2009) 



or 



that their density distribution lias led to a decelerating, 



thus stable, bubble. A more recent study (Kuiper et al 



2012 1, in which both the frequency-dependent treatment 
and gray flux-limited diffusion are applied to the same 
initial conditions, seems to confirm that the instability is 
an artifact caused by the gray approximation of the di- 
rect stellar radiation. In any case, ionization would also 
fill the stabl e radiation-driven bubbles of |Kuiper et al.| 
(2011 20121 with ionized gas, so that stable radiation 
pressure- driven outflows cannot account for the observed 
high-mass molecular outflows either. 

Another source of momentum is lin e-driven winds of 
high-mass stars (e.g. Castor et al.|1975 ). Indeed, some of 
the observations show high-mass outflows that are poorly 
collimated and see m to be consistent with being wind- 
blown bubbles (e.g. Shepherd|2005 Beuther & Shepherd 
|2005[ ). However, typical massive outflows have momen- 
tum fluxes that are two orders of ma gnitude higher than 
those produced by stellar winds (e.g. Richer et al.| 2000). 
Hence, stellar winds cannot be the primary drivers of 
high-mass outflows already for reasons of momentum 
conservation, and they are certainly not responsible for 
the more collimated outflows found in earlier pha ses of 
high-mass s tar formation (e.g. 



Beuther et al . 

Recently, Peters et al. ( 2010ap have shown 



2002) 



;hat loniza- 

tion feedback can drive bipolar molecular outflows. Con- 
trary to classical understanding of H n region growth, 
these simulations show that the H n regions around mas- 
sive protostars do not expand monotonically as long as 



the high-mass star is still embedded in an accretion flow. 
Instead, the H n region flickers whenever the ionizing 
radiation gets shielded by dense filamen ts formed in the 



gravitationa lly unstable accretion flow ( Galvan-Madrid 
et al. 20111. This rapid fluctuation of the thermal gas 



pressure drives bipolar molecular outflows parallel to the 
net angular momentum vector: the pressure increase dur- 
ing the ionization process drives the material away from 
the density enhancement in the plane of rotation, and 
the subsequent recombination lets the blown-out mate- 
rial become a molecular outflow. 

3. SYNTHETIC ALMA OBSERVATIONS 

We use synthetic ALMA observations of CO t o char- 
acterize the io nization-driven outflows found by |Peters| 



etaL J 2010a). We describe the input m odels (Sec 



tion 



3T ) and the CO emission maps (Section 3.2). 



3.1. 



Input Models 

The simulations use an extended version of the 
adaptive-mesh num erical hydrodynami cs code FLASH 
( Fryxell et al.||2000 ) with sink particles ( |Federrath et al. 



2010[ ) to follow the collapse beyond the formation of the 



first objects, allowing us to describe the early evolution 
of the disk that builds up, and our improved hybrid char- 



acteristics raytra cing method ( Rijkhorst et al.[2006 Pe- 
ters et aL]|2010a l that propagates the ioniz ing and non- 



ionizing radiation on the grid. We refer to |Peters et al.| 
( |2010a[ ) for details on the numerical algorithm and ini- 
tial conditions. A thorough description of the spatial and 
temporal H II region evolution in these simulations can 



be f ound in [Peters et al.| d2010a| ) , |Peters et al.| (|2010c| 
and Galvan-Madrid et al. (2011 ). 

We analyze two simulations. Both simulations start 
from a 1000 M , spherical molecular cloud of diameter 
3.2 pc. The cloud is initially in solid body rotation with 
angular velocity to = 1.5 x 10~ 14 s -1 . During the col- 
lapse of this massive core, a rotationally flattened struc- 
ture forms that quickly becomes gravitationally unstable 
and fragments. In one simulation (Run B), we form a 
full stellar cluster in this rotating accretion flow, includ- 
ing radiative feedback by both ionizing and non-ionizing 
radiation. In another simulation (Run A), also including 
radiation feedback, we artificially prevent the formation 
of all stars after the first one and run the simulation with 
a single star only. In Run B, three massive stars with 
masses around 20 M form, each creating an H n region 
around itself, along with numerous lower mass stars. In 
contrast, the star in Run A reaches more than 70 M©. A 
thorough discussion of the effects of radiative feedback 
and the suppression o f fragmentation can be found in 
Peters et alT] ( |2010a|b| ) . 

The massive stars in both simulations drive bipo- 
lar outflows with their ionization feedback. The ba- 
sic mechanism, as already mentioned in Section |2.2[ is 
the repeated shielding of the ionizing radiation by dense 
streams of gas in the accretion flow around the massive 
star. The ionizing radiation accelerates the molecular gas 
around the H II region down the steepest density gradi- 
ent, perpendicular to the rotational flattening. When the 
ionizing radiation gets shielded by the filamentary accre- 
tion flow, the ionized gas can recombine on timescales as 
short as ~ 100 yr, and then the outflow becomes molec- 
ular, given that the formation timescale of molecular hy- 
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drogen in these den se outflows is less than a year (Hol- 
lenbach et al.|1971 |. The molecular outflow continues to 
expand by momentum conservation. Since the outflow 
is not driven anymore without the thermal pressure of 
the ionized gas, the outflow can escape from the gravita- 
tional potential well of the stellar cluster only if its initial 
velocity permits. After the loss of thermal support, the 
outflow decelerates and eventually falls back onto the 
rotationally flattened structure b y its gravitati o nal at- 
traction. The online material of |Peters et al. (2010a I 
contains movies of density slices showing this process tor 
both runs. 

Generally, the outflows in Run A are much more pow- 
erful than the outflows in Run B. The reasons are the 
stronger accretion flow in which the massive star is 
embedded and the larger ionizing luminosity. Though 
Run A is not a realistic model of massive star forma- 
tion globally, the driving of the molecular outflow by 
the interaction of the ionizing radiation with the accre- 
tion flow is modeled correctly. These stronger accretion 
flows could be realized in simulations with multiple stars 
by, for example, starting from a more massive core, us- 
ing a more centrally concentrated density profile (e.g. 
Girichidis et al. 2012), or reducing the initi al rotational 
motion of the core Je.g. ISeifried et alp)!!] ). 

The online material for this article contains animations 
of the ionization-driven outflows for both simulations. 
The videos show volume renderings of gas density (left) 
and thermal pressure (right). The direct comparison of 
the two images allows the identification of the H n re- 
gions and illustrates how the thermal pressure of the ion- 
ized gas drives the molecular fountains. For comparison, 
we also show an anim ation of a third simulation, Run E 
(see Peters et al. 2011 ), in which an additional large-scale 
outflow is driven by magnetic fields and interferes with 
the ionization-driven outflows. Though the interplay be- 
tween the different types of outflows is highly complex, 
the morphological differences between the magnetically- 
driven outflows and the ionization-driven outflows be- 
comes very clear. 

3.2. CO Emission Maps 

We have previously used the three-dimensional 
adaptive-mesh radiative transfer code RADMC-3D^] to 
generate maps o f free-free and dust emission from our 
simulation data (|Peters et al.|[2010c[ ). Here we use the 
capability of RADMC-3D to model molecular line emis- 
sion (e.g. Shetty et al.||2011 ). In order to trace the bulk 
dynamics of the gas, we have chosen to model the CO 
J = 2 — 1 transition at 230.538 GHz. We have simulated 
observing the outflows at a distance of 1.3 k pc, the same 



distance as G5. 89-0. 39 ( |Motogi et aLpOllh . This source 
is believed to be powered by a 25 Mq star, and is so 
directly comparable with one of the stars in Run B. 

CO is an excellent tracer of the bulk gas dynamcis in 
star-forming regions because of its relatively low critical 
density, and relatively high abundance (10~ 4 with re- 

1994] ). It 
lerefore a 



spect to molecular hydrogen, Wilson & Rood 
is the most abundant molecule after H2 and is t 
good tracer of large-scale structures. Its low critical den- 
sity means that it is detected even in the tenuous high- 

8 http:/ /www. ita.uni-heidelberg.de/~dullemond/software/radmc- 
3d/ 



velocity gas in outflows. SiO, another excellent outflow 
tracer, was not used because it emits on much smaller 
scales and its emission is directly tied to the type of 
shocks produced in the outflow (i.e. Si is liberated from 
dust grains by sputtering during the extended passage 
through C-shocks me diated by ion-neutral drag, but not 
by classical J-shock s ( Schilke et al. 1997 Gusdorf et al. 
2008l|Guillet et al.||20ll| ). 



For simplicity, we assume local thermodynamic equi- 
librium (LTE) conditions. We use the CO Einstein coef- 
ficients from the Le iden Atomic and Molecular Database 
( Scho ier et al.|20 05 ) . These 3D data cubes were then con- 
verted into skymaps, and observations with ALMA were 
simulated usin g the "simobserve" an d "simanalyze" tasks 
in CASA 3.4 ( |McMullin et al.||2007[ ). The native resolu- 
tion of the original datacube was 49 AU. Assuming the 
source was at the distance of 1.3 kpc, this gives a native 
resolution of 0.038" . This is much higher than the final 
resolution achieved in our ALMA simulations. We sim- 
ulated observing with ALMA in a moderately high res- 
olution mode ("alma.outlO.cfg"), achieving a resolution 
of 0.5" at 230 GHz with a channel spacing of 1 km s . 
The rms sensitivity for the ALMA observations should 
be 1.2 mJy/beam for a 0.5" beam. Noise was added in 
the UV plane to simulate 4 th quartile observing condi- 
tions. We also simulated observing with CARMA in its C 
configuration. This resulted in a final image resolution of 
1.3" x 1.0". The CASA simulator assumed a system tem- 
perature of 200 K for the CARMA observations. For all 
simulations, we assumed a total integration time of four 
hours. For the ALMA mosaic, we required six pointings, 
and with CARMA, two. 

Cleaning was done non-interactively using natural 
weighting and a 3 a cutoff. Zcroth and first moment 
maps were created using line wing emission, excluding 
emission from within 3 km s _1 of the source rest veloc- 
ity. The zeroth moment (integrated intensity) maps were 
used to compute th e column density follow ing the pro- 



cedure described in Klaassen et al 
temperature of T 



2011), assuming a 
1U0K. Because the size of the emit- 



ting region can be determined from the images, the num- 
ber of molecules was calculated from the column density 
assuming a CO abundance of 10 -4 with respect to H2. 
From this, the gas mass was determined (see the first 
column of Table Rl). From the first moment (intensity 
weighted velocity) maps, characteristic velocities of the 
outflowing gas were derived. These values are presented 
in the second column of Table [TJ The uncertainties in 
Table [l] reflect the noise output from the ALMA obser- 
vation simulator which were propagated (in quadrature) 
through the equations used to derive the outflow prop- 
erties. 

Outflow momentum was calculated from P — mv, and 
energy from E = 0.5mv 2 . Outflow luminosities and mass 
loss rates were determined by dividing the outflow ener- 
gies and masses by the kinematic age of the outflows 
(400 yr) . More details on this type of calculation can be 
found in Klaassen et al. (2011 1. For each of our modelled 
outflows, we find that many outflow properties as mea- 
sured by ALMA are about a factor of two higher than 
those observed with CARMA. 

Figures [l] and [2] show three color images of the CO 
emission for snapshots of Run A and Run B as would be 
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TABLE 1 

Outflow parameters derived from CARMA and ALMA simulations 



CARMA 




M 


V 


P 


E 


L 


M 


R 






(M e ) 


(km s- 1 ) 


(Me km s _1 ) 


(10 44 erg) 


(Le) 


(10- 3 Me yr- 1 ) 


(AU) 


Run A 


blue 


0.059±0.001 


3.816±0.126 


0.224±0.010 


0.085±0.007 


0.175±0.048 


0.147±0.031 


4100 




red 


0.066±0.001 


3.391±0.150 


0.223±0.012 


0.075±0.007 


0.155±0.046 


0.165±0.035 


4100 


Run B (left) 


blue 


0.013±0.000 


3.634±0.199 


0.048±0.004 


0.017±0.002 


0.036±0.012 


0.033±0.007 


3300 




red 


0.060±0.000 


3.318±0.199 


0.198±0.013 


0.065±0.008 


0.134±0.044 


0.149±0.031 


2100 


Run B (right) 


blue 


0.058±0.000 


3.549±0.140 


0.206±0.009 


0.073±0.006 


0.150±0.043 


0.145±0.030 


5000 




red 


0.034±0.000 


2.453±0.000 


0.082±0.001 


0.020±0.000 


0.041±0.009 


0.084±0.018 


4100 


ALMA 




M 


V 


P 


E 


L 


M 


R 






(M G ) 


(km s" 1 ) 


(Me km s _1 ) 


(10 44 erg) 


(Le) 


(10" 3 M e yr -1 ) 


(AU) 


Run A 


blue 


0.082±0.000 


3.888±0.099 


0.320±0.009 


0.124±0.007 


0.255±0.065 


0.206±0.042 


4100 




red 


O.lOliO.OOO 


3.297±0.154 


0.333±0.016 


0.109±0.010 


0.225±0.066 


0.252±0.051 


4100 


Run B (left) 


blue 


0.050±0.000 


3.446±0.000 


0.171±0.001 


0.059±0.000 


0.121±0.025 


0.124±0.025 


3300 




red 


0.141±0.000 


2.757±0.184 


0.388±0.026 


0.106±0.014 


0.219±0.073 


0.352±0.071 


2100 


Run B (right) 


blue 


0.060±0.000 


3.550±0.143 


0.213±0.009 


0.075±0.006 


0.155±0.044 


0.150±0.030 


5000 




red 


0.044±0.000 


2.370±0.000 


0.104±0.000 


0.024±0.000 


0.050±0.010 


0.109±0.022 


4100 



Note. — The table shows the mass M, velocity V, momentum P, kinetic energy E, mechanical luminosity L, mass-loss rate M, and 
size R of the observed outflows for two snapshots in Run A and Run B. The kinematic age of all outflows is 400 yr. 



observed with CARMA and ALMA, respectively, after 
being processed through CASA. The images were ob- 
served at an inclination of 30° with respect to t he disk 
normal direction , which is a typical orientation (Cabrit 
& Bertout 1986). The zeroth moment maps of the blue 
and red shitted emission were cut at 3 a for this plot 
(as shown in blue and red, respectively). The white ar- 
eas correspond to the brightest regions in the map, the 
position of the high-mass protostar. The two outflows 
visible in the map for Run B are driven by two sepa- 
rated massive stars, one of which still has a small H 11 
region (right), while the H 11 region around the other 
star has totally recombined (left). As the stars move 
through the flattened structure, their ionizing radiation 
continuously drives material out perpendicular to it, but 
after the star has moved and the gas has recombined, 
the gas falls back. This process creates a large-scale 
molecular fountain, but its velocity is smaller than our 
low-velocity cut-off, and the interferometer additionally 
filters out some of its large-scale structure. Hence, the 
map only shows the smaller, more powerful outflows that 
are currently driven by the ionization feedback from the 
massive stars, and not the large-scale bubble that is also 
present. 

4. DISCUSSION 
4.1. Comparison to exact solution 

We can calculate the mass, momentum and kinetic 
energy in the outflows exactly. To this end, we use 
RADMC-3D to sum over the mass, momentum and 
energy in zones with velocity high enough to be con- 
tained in the outflow as viewed along the line of sight 
of the simulated observation. We only take cells with 
3km/s< v < 7km/s into account, which corresponds to 
the selection criterion used with the synthetic maps. The 
results of these direct measurements are listed in Tabled 
The directly measured values are comparable with the 
numbers inferred from the synthetic observations, but 
generally higher. The reason could be that the synthetic 
observations do not recover all the flux from the outflows. 

The morphology of the outflow for Run A as pictured 



in Figure [2] suggest that the outflow is pointing mostly 
along the line of sight to the observer, which can be con- 
firmed by examining the raw simulation data. This find- 
ing is also supported by the intensity-weighted velocities 
(characteristic velocities, listed in Table ll| of the gas in 
this outflow, which are slightly higher than those for the 
outflows in Run B. The morphology of the two outflows in 
the figure for Run B (along with their lower characteris- 
tic velocities) suggest they are more likely oriented closer 
to the plane of the sky, which is again confirmed by the 
raw data. However, the ionizing luminosity that powered 
these outflows before the H 11 region recombined is also 
vastly different: the star in Run A has 57 M Q , while the 
two massive stars in Run B have masses around 20 A/ Q 
each. 

Our analysis shows that the flickering of the ioniza- 
tion feedback can produce molecular outflows that are 
morphologically similar to but smaller than the observed 
ones. They appear to be very young compared to the 
age of the star because they are not driven continuously. 
Because of their short kinematic age, the momentum, en- 
ergy and luminosity of the outflows are relatively small 
compared to observed outflows around high-mass stars. 
We predict that these small, high-momentum, but low- 
velocity outflows will be found by upcoming ALMA ob- 
servations that probe the driving sources of high-mass 
outflows, if they can be disentangled from the larger high- 
velocity outflows. 

4.2. Comparison to observations 

The outflows in Figure [2] are smaller than those gen- 
erally seen in high-mass star forming regions. However 
their kinematic ages (~ 400 yr) suggest they are also 
much younger than those which have been observed to 
date. These two properties are presumably correlated. 
Although the calculated outflow masses and kinematic 
properties (momentum, energy and mechanical luminos- 
ity) are consistent with the range of these properties 
found for high-mass star-forming regions (e.g. Wu et al. 
2004J, they are located near their lower boundaries 



We find the outflow velocities in the simulations to 
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Fig. 1. — Three color images of synthetic CARMA observations of snapshots from Run A and Run B processed with CASA. The resolution 
of these observations is 1.3" X 1.0", and the two point mosaic was simulated to have been observed for four hours on source. Noise was 
added in the UV plane to simulate a system temperature of 200 K, and the source was assumed to be at the distance of G5.89 (1.3 kpc). 
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Run B (Left) 
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Fig. 2. — Three color images of synthetic ALMA observations of snapshots from Run A and Run B processed with CASA. The resolution 
of these observations is 0.5", and the six point mosaic was simulated to have been observed for four hours on source. Noise was added in 
the UV plane to simulate fourth quartile observing conditions, and the source was assumed to be at the distance of G5.89 (1.3 kpc). 



be an order of magnitude smaller than those in G5.89 
even though previous observati ons of this source (e.g. 
Puga et al.||2006j |Su et al.|[2012| were at similar angular 



resolutions. 

The simulated mass infall rates are high compared to 
observations if we assume all of the infalling mass makes 
it onto the central forming star (Mi n ~ 10 -3 Mq yr _1 
suggests a 10 Mq star forms in 10 yr), however much of 
the infalling material will not be accreted by the central 
star(s). Instead, secondary stars form in the gravita- 
tionally unstable accretion flow around the massive star 
and intercept gas that would otherwise fall inwar ds in a 
process called fragmentation-induce d starvation ( Peters 
et~al]|2010a|b| |Girichidis et al.|[20l2] ). 



5. CONCLUSIONS 

Our results show that ionization feedback is able to 
drive molecular outflows around high-mass stars that 
have properties consistent with some observed ones. 
However, the measured properties are rather on the lower 
boundary of the typical values of high-mass outflows, 
which makes it unlikely that ionization feedback is the 
main driver of high-mass outflows. Additional feedback 
mechanisms like radiation pressure and stellar winds or 



TABLE 2 

Outflow parameters derived from direct measurements 







M 


P 


E 






(Mq) 


(Mq km s _1 ) 


(10 44 erg) 


Run A 


blue 


0.67 


2.43 


0.91 




red 


1.49 


5.84 


2.44 


Run B (left) 


blue 


0.04 


0.18 


0.08 




red 


0.23 


0.90 


0.38 


Run B (right) 


blue 


0.13 


0.53 


0.23 




red 


0.04 


0.17 


0.07 



Note. — Outflow mass M, momentum P and kinetic energy E 
as directly measured from the simulation data. 



magnetic fields might help, but as discussed in the Intro- 
duction, there is no strong evidence that this is indeed 
the case. 

We note, however, that lower-mass companions nat- 
urally form with high-mass stars in their gravitation- 



ally unstable acc retion flow (Peters et al. |2010a|b 2011 
Girichidis et al. 2012), in agreement with the observa- 



tion that massive stars are typically found as mem bers 
of higher-order multiple stellar systems and clusters ( Ho 
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high-mass outflows. Their contribution is orders of mag- 
nitude larger than the momentum in the ionization- 
driven outflows. 

Further simulations that simultaneously include ion- 
ization feedback from massive stars as well as the me- 
chanical feedback from their low- and intermediate-mass 
companions are required to thoroughly test this proposal. 



0.60 0.65 0.70 

i(Myr) 



0.75 



Fig. 3. — Momentum input P according to the outflow model 
from the stars in the cluster as function of time t. The joint mo- 
mentum injected by the low- and intermediate-mass stars is well 
within the regime of high-mass outflows. 



& Haschick 1981 Zinnecker & Yorke 20071. We pro 



pose that it is the joint outflow activity of the low- and 
intermediate-mass companions that we are seeing. 

We tested whether the low- and intermediate-mass 
companions in the accretion flow around the massive 
stars in our simulations at least in principle are able to 
reproduce the characteristics of high-mass outflows with 
a simple model. We assume that a factor of / = 0.1 of 
the accreted material onto the stars is ejected again as 
an o utflow with an initial velocity of v W md = 150kms _1 
(e.g. Bontemps et al.|1996 ) . The momentum injected into 
the surrounding gas by a single star per unit time is then 
P — f v w i n dM acc for an accretion rate M acc . Figure 
shows the integrated total momentum P as function 
time for all stars in the cluster. 

The figure shows that even with conservative assump- 
tions the low- and intermediate-mass stars alone can eas- 
ily produce outflow momenta in the regime of observed 
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